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Abstract 

We suggest the application of nitronylnitroxide substituted with methyl group and 2,2,6,6- 
tetramethylpiperidin organic radicals as 1/2-spin qubits for self-assembled monolayer quantum 
devices. We show that the oscillating cantilever driven adiabatic reversals technique can provide 
the read-out of the spin states. We compute components of the g'-tensor and dipole-dipole inter- 
action tensor for these radicals. We show that the delocalization of the spin in the radical may 
significantly influence the dipole-dipole interaction between the spins. 
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1. INTRODUCTION 



The self assembled monolayer (SAM) molecular systems with no or small amount of de- 
fects are the promising candidates for many electronic devices [ij. SAM systems containing 
radicals with unpaired spin 1/2 can be used in quantum logic devices ^. It was suggested 
using micro- wires, which produce a magnetic field gradient, and rf pulses, which induce 
Rabi transitions, to manipulate with the radical spins and create the quantum entangle- 
ment between them. Later the conditions for the molecules appropriate for the quantum 
logic devices have been formulated, and 1,3-diketone radicals have been suggested and theo- 
retically analyzed from the point of view of their application to the quantum logical devices 

This paper pursues two objectives. First, we suggest two stable organic radicals, 
nitronylnitroxide substituted with methyl group (NITRO) and 2,2,6,6-tetramethylpiperidin 
(TEMPO), for SAM quantum devices. In order to stimulate the experiments we compute 
the (yf-tensor and zz-component of the dipole-dipole interaction tensor. Second, we sug- 
gest the novel oscillating cantilever driven adiabatic reversals (OSCAR) technique as the 
tool for read-out of the spin states. The OSCAR technique has been recently successfully 
implemented for the single-spin detection |^. In section II we describe the effective spin 
Hamiltonian for the two interacting radicals, in section III we consider the method of cre- 
ation and detection entanglement between the two spins, in section III we describe the 
methods of computation of the (7-tensor and the dipole-dipole interaction tensors, and in 
section IV discuss the results of our quantum chemical computations. 



2. SPIN HAMILTONIAN FOR TWO QUBITS 



Let us discuss the simplest possible quantum computer element - the two qubit system. 
In SAM the two organic radicals are embedded in a monolayer. Electron spins S = 1/2 
localized on the radicals represent the two qubit system. We use geometrical arrangement 
the same as in the Ref. 1.2,1 , where external magnetic field Bq is oriented along the 2;-axis and 
radicals are separated by the distance a along the y-axis (see Fig. 1). In addition to constant 
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external magnetic field Bq a gradient of external magnetic field dBz/dy along the y-axis is 
applied on the SAM. If radicals in SAM are oriented in such a way that the principle z-axis 
of the electronic g-tensor corresponds to the orientation of the external magnetic field Bq, 
then the effective Hamiltonian for the spin group is 

dB 

Hspin = fiB9zzBoSz{2) + HB9zziBo - -^a)S^{l) - D^^S^{l)S^{2) , (1) 

where ^^(l), 5*2(2) are respectively radical one and radical two spin projection operators, fiB 
- the Bohr magneton, Qzz - the component of the electronic g-tensor of a radical (we assume 
that the radicals are the same in the spin group) and D^z is the dipole-dipole interaction 
tensor component. The Hamiltonian (0) can be treated in electron spin functions basis 
WiO'2) (o"i, o"2 = or 1, where stands for the spin ground state and 1 for the spin excited 
state). As in our system the magnetic field gradient along the y-axis exists, both radicals 
in the spin group (first and second) have distinguished spin transition frequencies, which 
depend on the local magnetic field strength at the radical position and additionally shifted 
by the dipole-dipole interaction between the radicals. 

Under these conditions, energy levels of the system are defined by the following expres- 
sions: 

|00) : E = -^fiB9zzBo - ^f^BQzziBo + Ga) + ^D, (2) 

|01) : E = ^nsgzzBo - ^HBgzziBo + Ga) - ^D, 

|10) : E = -^fiBgzzBo + ^fiBgzziBo + Ga) - ^D, 

: E = ^fiBgzzBo + ^fiBgzz{Bo + Ga) + ^D, 

G = \dBz/dy\, D = \Dzzl 

{Dzz and dB^/dy are negative), the external magnetic field on the second spin is Bq and on 
the first spin is Bq + Ga). The energy levels diagram of two qubit system with four possible 
transitions is given in Fig. 2. The transition frequencies are determined by the following 
expressions {h = 1) : 

|00) - |01) : = fiBgzzBQ - ^D, (3) 
|00) - |10) : uj',= iXBgzziBQ + Ga) - ^D, 
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|10)-|11): cul = iXBgzzBo + lD, 

|01) - 111) : ul = f^BQzziBo + Go) + ]^D. 

Here cuf means the transition frequency for the spin "i" {i = 1,2) if the neighboring spin 
is in the state \k) {k = 0,1). Given formulas for radical spin energy levels involve two 
molecular parameters: electronic g-tensor component g^z (property of radical itself) and 
dipole-dipole coupling parameter D, which is the property of the two radical system. Both 
these properties can be efficiently evaluated using quantum chemistry methods and influence 
of various parameters (radicals geometries, distance between radicals in spin group, radical 
relative orientation to each other, etc.) on these properties can be investigated. 



3. ENTANGLED SPIN STATES 



In this section, we consider the creation and detection of the entangled state 

-^[|OO)+exp(i0)|ll)], 

where is a nonsignificant phase factor. 

First, one applies a 7r/2-pulse with the frequency cul (we assume that initially both spins 
are in their ground state |00), as it is shown in Fig. 1). This pulse drives the spin system 
into the state 

-^[|0) + exp(i</.0|l)]|0). 

Now one applies a 7r-pulse with the frequency luI- This pulse drives the second spin if the 
first spin is in the state 1 1) . Thus, the new state of the system is 

-^[|OO)+exp(z02)|ll)]. 
This is a desired entangled state of the two-qubit system. 

To detect this state one uses the OSCAR technique. In the OSCAR technique the 
cantilever tip with the ferromagnetic particle oscillates along the line which is parallel to the 
sample surface. The distance between the ferromagnetic particle and the selected radical 
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spin changes periodically. When the cantilever tip is in its end point one starts to apply 
the rf field. The magnetic field on the spin changes its magnitude with the cantilever 
period. Suppose that the equilibrium position of the cantilever corresponds to the resonant 
condition for the spin. Then, in the rotating system of coordinates the effective field on the 
spin changes its direction from +z to -z in the x-z plane (we assume that the rotating r/ field 
points in the positive x-direction of the rotating system). The condition of the adiabatic 
reversals is 

« flBQeBt, (4) 



' dt 

where B^f is the effective magnetic field in the rotating frame, Bi is the r/ field amplitude, 
and Qe is the free electron g-factor. If the condition of the adiabatic reversals is satisfied the 
spin follows the effective field. 

The back action of the spin on the cantilever tip causes the frequency shift of the can- 
tilever vibrations, which can be measured with high accuracy. The cantilever frequency shift 
can take two values ±|5ct;c| depending on the spin direction relatively the effective field. Mea- 
suring the sign of the frequency shift one can find the initial direction of the spin relatively 
to the effective magnetic field. We assume that initially (when we start application of the rf 
field) the effective magnetic field has the same direction as the external magnetic field. In 
this case the ground state of the resonant spin corresponds to the negative frequency shift 
of the cantilever vibrations Q . 

How to verify the entangled state of the two-spin system? The entangled state manifests 
itself in the two outcomes of the measurement: |00) or |11). Let, for example, we apply 
the r/ field with the frequency u = 002- In general, we have three possible outcomes for the 
measurement of the cantilever frequency shift: 1) The cantilever frequency shift is negative. 
It means that the second (right) spin is in the ground state, and the first (left) spin is 
also in the ground state. Thus, our system collapsed to the state |00) (in the system of 
coordinates connected to the effective field). 2) The frequency shift is zero. It means that 
our system initially collapsed to the states |11) or |10) (or their superposition), which are 
insensitive to the frequency 00 = cug- In this case, after the integer number of the cantilever 
cycles we may change the r/ field frequency from u = to, for example, ut = u^. Now we 
have two possible outcomes: 2a) For the state |10) the left resonant spin will provide the 
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positive cantilever frequency shift; 2b) For the state |11) there is no resonant spin, and the 
frequency shift is zero. 3) The frequency shift is positive. It means the second resonant spin 
is in the excited state |1), and the first spin is in the ground state |0), i.e. our system has 
collapsed to |01). Thus, the repeated outcomes 1) or 2b) (with equal probability) correspond 
to the initial entangled state, which collapses to the states |00) or |11). The outcomes 2a) 
or 3) will show that the system is not in the expected entangled state. One can see that 
the OSCAR technique provides the simple verification of the quantum entanglement. Note 
that the verification experiment must be conducted for the time interval smaller than the 
characteristic time between the spin quantum jumps. 



4. COMPUTATIONAL METHODS 



A. Electronic (7-tensor 

According to Spin Hamiltonian (see Eq. ^ electronic (7-tensor of a radical is defined as 
the second derivative of the molecular electronic energy E 

g = T^^T^ ■ (5) 

Hb oBdS B=o,s=o 

For molecules described by a Breit-Pauli Hamiltonian, where the spin and the magnetic field 
as well as all other relativistic corrections are treated in perturbation theory framework, the 
molecular electronic (^-tensor can be evaluated using expression (up to second order in the 
fine structure constant a) 0|: 

g = gel + AgRMC + AgGC(le) + AgGC(2c) + Agoz/SO(le) + Agoz/SO(2e) • (6) 

In this equation, the first terms is the free electron g-factor, which comes from the electronic 
Zeeman operator in the Breit-Pauli Hamiltonian, with the radiative corrections from quan- 
tum electrodynamics introduced into the Hamiltonian empirically. The next three terms 
originate from first order perturbation theory applied to the Breit-Pauli Hamiltonian and 
are the mass-velocity and the one- and two electron corrections to the electronic Zeeman 
effect. The last two terms are, respectively, the one- and two-electron corrections contribut- 
ing to the electronic g-tensor to second order in perturbation theory as cross terms between 
the spin-orbit operators and the orbital Zeeman operator. All terms except the free electron 



g-factor value, Qe, in Eq. © contribute to the electronic g-tensor shift Ag, which accounts 
for the influence of the local electronic environment in the molecule on the unpaired elec- 
trons. From the contributions to the g-tensor shift given in Eq. (jH)), the first three terms 
are evaluated straightforwardly according to Eq. (0) from expectation values of the corre- 
sponding Breit-Pauli Hamiltonian operators. The last two terms, namely spin-orbit (SO) 
corrections, are considerably more involving computationally, as they are defined in terms 
of second-order perturbation theory, i.e. their calculation formally require a knowledge of 
all relevant excited states energies and wave functions. The evaluation of these terms us- 
ing density functional theory (DFT) methods are most often done usins various kinds of 



approximate sum-over-states approaches or response theory methods 



In the following part of this section we briefly describe electronic (7-tensor evaluation 



approach implemented in quantum chemistry code DALTON 10], which is based on a 
restricted DFT linear response theory ^. In this approach from the first-order perturbation 
theory contributions to the electronic g-tensor shift, we only include those terms that involve 
one-electron operators as DFT in principle can 
Cartesian ah components of these contributions to Ag tensor are 

Ag^W = -y(0|EP^^.|0)5'^', (7) 
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where Pi is the canonical linear momentum of electron Sz^i the z-component of the spin 
operator of electron z, rj^r and are the position vectors of electron i relative to nucleus 
and the magnetic gauge origin O, respectively. In the above equations \l/o is chosen to 
be the ground-state wave function with maximum spin projection S = Ms- Neglecting of 
the two-electron gauge correction, AgGc(2e), the contribution to the gf-tensor shift in DFT 
calculations does not influence the accuracy of the g-tensor evaluation as this term only gives 
a c„„ec«o„ f..o. a two-eleCon sc.eeuiug of the Ag.c,.. a„d ^considering the sn,a..ness of 



the one-electron gauge correction term itself is justified p, |7|, |9|. The major contributions 
to Ag tensor arising from second-order perturbation theory, so-called one- and two-electron 
SO corrections, are evaluated as linear response functions 

'^goz/so(ie) = ^((^io! -^so(io)))o, (8) 
7 



^g0Z/S0{2c) — ^((^io! -^SO(2o)))o, 

where the spectral representation of the hnear response function at zero frequency for two 
arbitrary operators Hi and H2 is given by 

{{Hi;H2))q= 2^ ^ . (9j 

m>0 " ^rn 

In Eq. (8) I^q is the cartesian a component of the angular momentum operator of electron 
and -f^so{ie) -^so(2e) Cartesian h components of the one- and two-electron SO 

operators. The Cartesian component h of the one-electron SO operator is defined as 



As mentioned above in DFT two-electron operators can not be evaluated properly and one 
need to introduce one or another approximation of the two-electron operators in order to 
perform calculations within limits of formalism. For g-tensor calculations performed in 
this work we selected to approximate two-electron SO operators by an Atomic Mean Field 
(AMFI) SO operator ^ as previous experience with the AMFI SO approximation in ah intio 
and DFT works devoted to electronic g-tensors calculations have been very encouraging and 
no significant problems with the accuracy of the AMFI SO approximation has been reported 

0. 

B. Dipole-dipole coupling 

Formally, in the case of two radicals with the spins S = 1/2, the interaction operator and 
corresponding D tensor may be written as 

rr 2 2 2 3(Si ■ r) (S2 " r) - r^Sj ■ S2 , . 

HDD = af^B9e > (11) 

Da = a fiB9e{ — :^ 

Dij = a'^fxlgji^^), {i,j = x,y,z), (i^j), 

where Si, S2 are respectively radical one and radical two effective spin operators, r - the S2 
position vector with respect to Si. In above definition of the dipole-dipole interaction tensor 
we assumed that each radical electronic g-tensor is isotropic and equal to the free electron g- 
factor. In this work we selected two methods for calculation of the dipole-dipole interaction 
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tensor required for determination of transition frequencies in the two qubit system (see Eqs. 
(3)): classical point dipole approach and "point dipole-spin density" approach. In case of 
classical point-dipole approximation delocalization of unpaired electrons is neglected and D 
tensor components are evaluated (in our two spin system geometry, r=(0,a,0)) as 

D,, = D,, = -a^^lgl^, (12) 

CL 

n o 2 2 2 l 

Uyy = 2a l^B9e-^^ 

One can expect this approach to be accurate for large separation of radicals, then the un- 
paired electrons delocalization region are very small comparing to the distance between the 
radicals. However, in case of moderate separation between radicals (0.5-1.5 nm) unpaired 
electrons delocalization in radicals can not be neglected and should be explicitly taken into 
account evaluating spin-dipole interaction tensor D. The simplest approach, which partially 
accounts for unpaired electrons distribution in radicals, is "point dipole-spin density" ap- 
proach in which one radical unpaired electron treated as delocalized and another radical 
electron magnetic moment represented as a point dipole. In this approximation D tensor 
components are evaluated as 

p{r,)±^-^dV, (13) 

Dij = a^B^e / Pi^e)^^dV, = x,y,z), {i ^ j), 

where p(re) is the electron spin density distribution in radical and r is the electron position 
with respect to the magnetic dipole position fm, r — Vg — tm- The described approach only 
partially takes into account delocalization of the unpaired electrons in D tensor calculations, 
but in order to obtain qualitative picture of the electron delocalization influence on the 
dipole-dipole interaction of the two radicals is sufficient. In this work we use this approach 
to determine limits of the classical point dipole approach for evaluation of the dipole-dipole 
interaction tensor D in two radicals system. 
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5. COMPUTATIONAL DETAILS 



As the promising candidates for qubits in SAM of organic radicals we selected two well- 
known stable organic radicals, namely methyl group substituted nitronyl nitroxide (NITRO) 
and 2,2,6,6-tetramethylpiperidin (TEMPO) (see Fig. 3), which form stable organic crystals 
and already have applications in material sciences as well as biochemistry. Both chosen 
radicals possess one unpaired electron and therefore each radical can be recognized as a single 
qubit. The geometry of NITRO and TEMPO radicals used in calculations of electronic g- 



tensors and of dipole-dipole interaction tensor D have been obtained by per: 
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optimization of single radical in 6-311G(d,p) basis set 
All geometry optimizations have been carried out in Gaussian-98 program 
optimizing geometries of NITRO and TEMPO radicals we performed geometry optimization 
of the NITRO and TEMPO radicals derivatives with (CH2)n tails (see Fig. 4). The (CH2)n 
groups were added to radicals in order to simulate conventional structure of the compounds 
used in formation of SAM, which usually feature long (CH2) groups or similar tails. Obtained 
structures of NITRO and TEMPO radicals with (CH2)n tails have been used to investigate 
influence of the tails on the properties of the NITRO and TEMPO derivatives compared to 
free NITRO and TEMPO radicals and allowed us to estimate feasibility of the NITRO and 
TEMPO radicals usage as basic building blocks of the spin arrays in SAM. 

Calculations of electronic g-tensors for NITRO and TEMPO radicals as well as their 
derivatives with (CII2) tails have been carried out using BP86 exchange-correlation func- 



tional [13|, ll7[, which gives accurate results for organic radicals. In all calculations of g- 
tensors we employed IGLO-II basis set ^| especially designed for evaluation of magnetic 
properties. In order to estimate influence of the (CH2)n (n=l,2) tails, which usually added 
to radicals in order to enable formation of SAM, on the electronic g-tensors of the radicals 
we carried out electronic g-tensors calculations for single NITRO and TEMPO radicals (see 
Fig. 3) as their derivatives with (CH2)n tails (see Fig. 4). 

Calculations of the dipole-dipole interaction tensor between two NITRO or two TEMPO 
radicals have been performed using previously described point dipole and "point dipole- 
spin density" approaches varying the distance between radicals, a, from 1 nm to 2 nm. 
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Single radicals optimized geometries have been employed in calculations using "point dipole- 
spin density" approach. All calculations have been performed at the B3LYP level using 

n 

Duning's double zeta basis set [191], which allow adequate description of the electron density 
distribution in investigated radicals. In calculations of the dipole-dipole interaction tensor 
we limited ourselves only by investigation of the D^z component dependence on the distance 
between radicals a and contrary to investigation of the electronic g-tensors do not carried 
out investigation of the (CH2)n tails influence on dipole-dipole interaction tensor D, as our 
g-tensor calculations indicated negligible influence of the (CH2)n tails on unpaired electron 
density distribution in both NITRO and TEMPO radicals. 



6. RESULTS AND DISCUSSION 



A. Electronic g-tensor 

Electronic g-tensor calculations results for single NITRO and TEMPO molecules as well as 
their derivatives with (CIl2)n n=l,2 tails are tabulated in Table l.nThe results of electronic 
g-tensor calculations for NITRO compound with (CH2)n tail separated for two different 
conformations of this compound, which differ only by orientation of the tail with one (CH2) 
unit. However, already for tail consisting of the two (CH2) units, there is no difference 
between NITRO A and NITRO B conformations due to the increased flexibility of tail, 
which leads geometry optimization procedure converges to same structure independently on 
the starting geometry. Here, we note even thought our calculations predict small radical 
tail rotation around C-C bonds it does not correspond to the "real" behavior of the (CH2) 
tail in SAM, as the motion of the tail is constrained by surrounding molecule tails in SAM. 
Now let us turn discussion from radicals geometrical structure features to their electronic 
g-tensors, which are one of the key quantities in our two qubit system Spin Hamiltonian. 

All electronic g-tensor components, presented in Table 1, B^re only slightly altered by 
addition of (CH2)n tail chain units, suggesting small distortion of unpaired electron density 
in radical by such chemical modification. The electronic g-tensor shifts almost converge, and 
extension of the chain further from two to three (CH2) units changes g-shift components 
in range 1-5 ppm. Therefore, we conclude that electronic g-tensors of the NITRO and 
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TEMPO radicals are only slightly affected by chemical addition of the (CH2)n tail, which 
is required for growth of the SAM, i.e. both radicals preserve their properties in SAM. 
Selected radicals have highly anisotropic electronic g-tensor (see Table 1 □ ) with large /^.g^^ 
and ^Qyy components (for orientation of the electronic g-tensor axes see Fig. 1). Only /S^g^z 
is small, and therefore along this axis we have molecular g-tensor component close to the free 
electron g-factor. The large anisotropy of the electronic g-tensor poses significant difficulties 
for the quantum computation prospective as in order to have well defined frequencies of the 
transitions it is essential to have fixed orientation of the electronic g-tensor principal axes 
with respect to the external magnetic field, i.e. radicals in SAM should be fixed in their 
position and rotations of the radicals with respect to their g-tensor principle axes must be 
restrained. 

B. Dipole-dipole interaction 

Dipole-dipole interaction tensor between two radicals is one of the key parameters in the 
spin Hamiltonian, which describes the two qubit system. We employed different approaches, 
namely classical point dipole approach and "point dipole-spin density" approach, for eval- 
uation of the D = iD^^l.The second approach as discussed previously partially accounts 
electron distribution in molecule and therefore results are directly dependent on the un- 
paired electron density localization in the molecule. Contrary, the first approach does not 
account for unpaired electron delocalization in radical and therefore gives the same results 
for all radicals. The dipole-dipole interaction dependence on the distance between radi- 
cals a is plotted in Fig. 5, where results for both classical point dipole approach (denoted 
classical) and "point dipole-spin density" approach (denoted Nitro and Tempo for NITRO 
and TEMPO radicals, respectively). Quick inspection of these plots indicates a substantial 
influence of the unpaired electron delocalization on the D value especially at the small a 
region. Therefore, unpaired electron delocalization can not be neglected in evaluation of 
the D tensor and previous data, which have been used in modeling of the qubits system, 
obtained with classical point dipole approach should be carefully reexamined for a range of 
the 10-15 A. Another implication of the non-negligible contribution from unpaired electron 
delocalization to the dipole-dipole interaction is that orientation of the radical can influence 
D tensor components values. Therefore, similarly to electronic g-tensor calculations the 
dipole-dipole interaction modeling results suggest that the rigid fixation of the radical in 
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SAM with specific orientation is one of the key requirements in production of SAM in order 
to make them useful in quantum computing. 

7. CONCLUSION 

We have suggested using NITRO and TEMPO radicals as spin qubits for a quantum 
logical device based on SAM systems. In order to stimulate experimental implementation of 
our idea we have computed the components of g-tensor and dipole-dipole interaction tensor 
for these radicals. We have shown that adding {CH2)n tail chain does not infiuence signifi- 
cantly the radical electron g-factor. Delocalization of the electron spin in radicals infiuences 
the dipole-dipole interaction between the radicals. We suggested a scheme for detection of 
the entanglement between the two radicals based on the novel OSCAR technique, which has 
been recently used for the single-spin detection. 
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TABLE I: Shift of the electronic g-tensor components on (CH2)n chain in nitronyl-nitroxide and 
TEMPO radicals. "'^ 





NITRO A 


NITRO B 


TEMPO 


Tail 


— CH2 (0112)2 


— CH2 (0112)2 


— CH2 (0112)2 




63 56 51 


63 59 51 


27 58 86 


Acj'22 


8922 8905 8882 


8922 8889 8838 


7289 7278 7245 


A533 


4206 4203 4206 


4206 4203 4206 


3945 3989 3791 



" Electronic g-tensor RDFT-LR calculations performed using BP86 exchange-correlation 
functional in Huz-11 basis set. Electronic g-tensor shifts are in ppm. 
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FIG. 1: Two spins S{1) and S{2) separated by the distance a. The external magnetic field B{y) 
is oriented in the positive 2;-direction and has a gradient in the y-direction. In the ground state, 
shown in the figure, the spins point in the negative z-direction, G = \dBz/dy\. 
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FIG. 2: The diagram of the energy levels for the two qubit system. 
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FIG. 3: NITRO and TEMPO radicals. 
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FIG. 4: NITRO and TEMPO radicals with CH2 tails (n = 1,2). 
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FIG. 5: Dependence of the parameter D on the distance between two radicals in the spin group (1 
- two interacting TEMPO radicals, 2 - two interacting NITRO radicals, 3 - two interacting point 
dipoles. 
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